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Fermilab experiment E835 has studied reactions pp — * tt°ty° , 7r ?i, -qn, -K°n' and r\r\ in the energy 
region of the Xco(l 3 Po) from 3340 MeV to 3470 MeV. Interference between resonant and continuum 
production is observed in the 7r°7r° and nn channels, and the product of the input and output 
branching fractions is measured. Limits on resonant production are set for the and ir°T]' channels. 
An indication of interference is observed in the nrj' channel. The technique for extracting resonance 
parameters in an environment dominated by continuum production is described. 

PACS numbers: 13.25.Gv;13.75.Cs;14.40.Gx 



I. INTRODUCTION 

E835 has measured the cross section of antiproton- 
proton annihilation into two pseudoscalar mesons (pp — > 
P1P2) m the energy region of the XcO- The final states 
studied are: 

P X P 2 = tt ^ (1), A (2), W (3), A' (4), and W (5). 

A total integrated luminosity of 33 pb -1 has been col- 
lected at 17 energy points from 3340 MeV to 3470 MeV. 
The mesons were detected through their decay into two 
photons. The cross sections reported here have been cor- 
rected for the respective branching ratios [1]: B(tt° — > 
77) = (98.798 ± 0.032)%, B{r) -> 77) = (39.43 ± 0.26)%, 
and B(t]' -> 77) = (2.12 ± 0.14)%. 

The analysis of reaction (1) has been published in letter 
format [2]. In the present work, processes (2)-(5) are 
reported for the first time and a more extensive discussion 
is given of process (1). More details on the analyses (1)- 
(3) can be found in a dissertation [3] . 

The expression for the angular distribution of pp — > 
P1P2, in the vicinity of the \c0 resonance, is given in 
Equations (6), (7) and (8). The subsequent isotropic 
decay of each meson into two photons need only to be 
considered for acceptance determination. 
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where 9* is the production polar angle of the two mesons 
in the center of mass (cm), while x is defined as 
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The term Ae l&A is generated from the pp state with he- 
licity A,; = A p — Xp = (hclicity-0), where X p and Ap are 
the proton and antiproton helicities, respectively. The 
helicity-0 coefficients and phases of the expansion are Cj 
and 5 j, respectively, while Pj(z) are Legendre polynomi- 
als. The term Be lSe is generated from the pp state with 
helicity |Ai| = 1 (helicity-1). The helicity-1 coefficients 

respectively, while P}{z) are 



and phases are Cj and Sj, 



Legendre associate functions Pj 1 with M = 1. The XcO 
(l 3 Po) resonance, parameterized by a Breit-Wigner am- 
plitude —Afj/(x + i), has been extracted from the J = 
term. Both nonresonant terms, Ae l&A and Be l&B , have a 
slowly varying (on the scale of the \c0 width, « 10 MeV) 
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FIG. 1: The first Legendre polynomials Pf =0 = Pj (left) 
and associated functions Pj^ =1 (right). 



dependence on x and an angular dependence on z de- 
scribed by the Legendre polynomials. The amplitudes 
— Aft/(x + i) and Ae lSA , from the same helicity-0 ini- 
tial state add coherently and interfere with each other. 
The amplitude Be lSn is added incoherently because it is 
generated from the helicity-1 initial state and does not 
interfere with the helicity-0 terms. 
Equation (6) can also be written as 

da A% , , sin 5 a — x cos 5 a o , 

— = -2~tt +A 2 + 2A R A ^— ± +B 2 . 11 

dz x z + 1 x z + 1 

S v ' 

interference— term 

Given the moderate energy dependence of Ae l5A and 
Be lSB , Equation (11) shows that at a fixed value of z, as 
E cm varies across the resonance, even a small resonant 
contribution can generate a significant interference signal 
superimposed on the continuum. As an example, for a 
resonant amplitude A R that is one tenth of the helicity- 
continuum amplitude A, the peak contribution to the 
cross section from the Breit-Wigner term A R /(x 2 + 1) 
is only 1% of A 2 , while the factor 2ArA of the inter- 
ference term of Equation (11) is 20% of A 2 . The factor 
(sin(>A — xcos5a)/(x 2 + 1) determines the shape of the 
interference pattern. A small resonance and a large con- 
tinuum is a typical experimental condition found when 
studying charmonium in the process pp — > hadrons. 

The term B 2 provides an additional contribution to 
the cross section, and if B 2 is too large compared to 
the interference term, it may mask the presence of the 
resonance. It is quite helpful that the associated Leg- 
endre functions with M = 1, which are the constituents 
of the amplitude Be %&B , share a common multiplicative 
factor z which causes them to vanish at z = (see Fig- 
ure 1). Since the polynomials Pj are either squared or 
multiplied by each other, B 2 is suppressed with respect 
to A 2 by z 2 at small z. To measure the value of the 
resonant amplitude An, it is critical to know the size of 
both the hclicity-0 continuum, with which the resonance 
is coherent, and the helicity-1 continuum, which is inco- 
herent with the resonance. The region at small z (i.e. 
cos 6** w 0) is therefore the natural place to concentrate 
for the determination of Ar. 



TABLE I: How the orbital angular momentum L and the spin 
S of the pp in the initial state combine to make different J PC 
values and the corresponding charmonium resonances. A / 
indicates whether a pseudoscalar-pseudoscalar (P1P2) state 
is accessible. Spectroscopic notation is given for pp ( +1 Lj) 
and cc (n 2S+1 Lj). 
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Table I shows how the orbital angular momentum 
(Lpp) and spin (Sp p ) combine to form different values of 
the total angular momentum, parity and charge conjuga- 
tion (J PC ). For each J PC , it shows which charmonium 
resonance is formed, and whether the J PC — even ++ 
pseudoscalar-pseudoscalar final states can be accessed. 
The spectroscopic notation ( 2S+1 Lj) is given for pp and 
cc states. Both states are fermion-antifcrmion, thus are 
described in the same way as far as L, S and J PC are 
concerned, while a radial quantum number (n) applies 
only to the bound cc system. The \c0 (J PC = ++ ) can 
be produced only from the spin-triplet L Pp = 1. Table I 
is sorted by increasing Lp p , which correlates with the pp 
impact parameter (b). Even for nonresonant reconfigura- 
tion of the pp system into two pseudoscalar mesons, small 
impact parameter is favored at small z since the valence 
quarks must either annihilate or suffer large momentum 
transfers. Excellent fits to the angular distributions of 
the two pseudoscalar mesons are obtained with a limited 
number of partial waves. 



II. DATA SELECTION 

A comprehensive description of the E835 apparatus 
and experimental technique is provided in [4]. The data 
selection and the determination of the acceptance and ef- 
ficiency for the processes (l)-(3) are extensively discussed 
in [3]; only a summary is given here. Processes (4) and (5) 
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are discussed in Section X. 

The stochastically-cooled p beam circulating in the an- 
tiproton accumulator intercepts a hydrogen gas-jet tar- 
get. The energy of the beam can be tuned to the energy 
of interest. Before the year 2000 run, the accumulator 
transition energy was raised to E cm of 3600 MeV. A tech- 
nique was developed to modify the accumulator lattice 
and lower the transition energy as the beam was decel- 
erated [5], thus allowing adequate margin between the 
operating energy and the transition energy. The E cm 
spectrum is approximately gaussian and is determined 
from measurements of the p-beam revolution frequency 
and orbit length. The precision of the measurement of 
the mean value of the E cm spectrum is about 100 keV. 
The r.m.s. spread of the E cm spectrum is a few hundred 
keV, much smaller than the width of the Xco resonance. 

The important detectors for this analysis are the cen- 
tral calorimeter (CCAL) which is used to measure the 
photon energy deposits, the system of scintillation coun- 
ters which vetoes events with charged particles and the 
luminosity monitor. The energy resolution of CCAL is 
a E /E ~ 6%/y/E(GeV) + 1.4%, while the polar and 
azimuthal angular resolutions are ag ~ 6 mrad and 
Ofp — 11 mrad, respectively. 

Online, two-body candidate events were selected by 
means of two independent triggers based on CCAL: 
the two-body trigger and the total-energy trigger. The 
two-body trigger accepted events with two large energy 
deposits in CCAL satisfying two-body kinematics; the 
total-energy trigger accepted events where at least 80% 
of the center-of-mass energy was deposited in CCAL. 
For 7r°7r° events, the efficiency of the two-body trigger 
is 99.9% - the small opening angles of the decay ir° — > 77 
keeps the tt°7t° — ► 27 + 27 events within the two-body 
trigger requirements - while the efficiency of the total- 
energy trigger is ~ 98.2%. The two-body trigger is some- 
what less efficient for 77 — > 77 decays and only the total- 
energy trigger was used for selecting ir°r] and 7777 events, 
with an average efficiency of ~ 98.5% for both channels. 
Both triggers were subjected to the charged-particle veto. 

Offline, events with four CCAL energy deposits greater 
than 100 MeV are retained if they passed a 5% probabil- 
ity cut on a 4-constraint (4C) fit to pp — > 47. There 
are three ways to combine the four photons into two 
pairs; the event topology was chosen as the combina- 
tion with the highest confidence level of a 6C fit to 
pp — > P1P2 — ► 47. However, given the limited open- 
ing angles of the decay 7r° (or 77) — > 77 at these energies 
and the large angle between the two mesons, it is virtu- 
ally impossible that more than one combination satisfy 
the 6C fits for more than one of the processes (l)-(3). 
Figure 2 is a LEGO plot of the two 2-photon invariant 
masses prior to any fitting. Evident are high peaks of 
7r°7r°, tt°t] and 7777 events, as well as those of ir a u> (where 
a photon in the oj — > ir°j — > 777 decay chain is not ob- 
served), n°ri' and 7777'. Marked "berms" are produced by 
ir°X events, where X can be one or more particles that 
partially escaped detection. Softer r\X berms are also 




FIG. 2: LEGO plot of the two 2-photon invariant masses 
m 7172 and m 7374 . Note the logarithmic scale of the vertical 
axis. 

noticeable. Prior to a 4C fit, the mass resolution for a 
7T°, an 77, and an 77' is cr~15 MeV, 25 MeV and 35 MeV, 
respectively, with a small dependence on z. 

A set of kinematic cuts is applied for each process with 
the values shown in Table II. These include cuts on the 
invariant masses of each pair of photons and on the co- 
linearity and coplanarity of the candidate mesons. The 
colinearity, AO, of the two mesons A and B, is given by 
^mcas ™ wnere $meas is the measured polar angle of 
meson A while 6^ lc is the value of the same quantity 
calculated assuming two-body kinematics using the mea- 
sured polar angle of meson B. In the case where the two 
mesons are the same, A is chosen to be the one in the 
forward direction. In the case where they are different, A 
is chosen to be the lighter of the two mesons. The copla- 
narity, A<f>, is defined as 180° - K loas - <^ eas |, where 
'/'meas an d ^mcas are tne measured azimuthal angles of 
the mesons. 

The selection efficiency (e) and detector acceptance (a) 
are determined by Monte Carlo (MC) simulation. MC 
events are generated with a uniform angular distribution 
and a matrix is built to record both the generated and 
reconstructed values of z. The angular distribution of the 
data is determined by using the inverse of this matrix to 
correct for effects due to the angular resolution of the 
detector. The product of the acceptance and efficiency 
(axe) as a function of the generated value of z for 7r°7r°, 
7r°?7 and 7777 is shown in Figure 3. 

A significant source of inefficiency is event pileup. It 
is particularly important to correct for this effect since it 
is rate dependent and the instantaneous luminosity var- 
ied from one energy point to another. To determine the 
pilcup effects, MC events were overlaid with the contents 
of random-gate events recorded throughout the data- 
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TABLE II: The kinematic cuts applied to each channel. The asymmetry of the cut on the invariant mass reconstructed of the 
candidate tt° compensates for a slight overestimate of the invariant mass reconstructed for 7r°'s when the energy deposits from 
the daughter photons overlap in the CCAL (see details in [4]). 
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FIG. 3: The product of detector acceptance and selection efficiency, axe, as a function of z for 7r ( tv" , it r) and rjri events. In 
the ir°ri case, where the two mesons are distinguishable, 2 is defined as cos 0*0 . 



taking, thus reproducing the conditions of each energy 
point. The event reconstruction was then performed on 
these hybrid events to determine the reconstruction effi- 
ciency. The rate-dependent losses vary from 14% to 23% 
with an average of 20%; these losses do not differ signifi- 
cantly among the analyzed reactions. As a function of z, 
the pileup correction is determined and applied. 



III. SIGNAL AND BACKGROUND 
SUBTRACTION FOR ttV , n°ri AND rjr) 
SAMPLES. 

Figure 4 shows the region of the 7r°7r° peak before ap- 
plying the mass cuts. A log-likelihood fit is performed 
to this data. The function used is the sum of a two 
gaussians, having the same mean, to describe the 7r°7r° 
peak, and two gaussian berms and a tilted plane to de- 
scribe the background (from events such as 7r°w — > 7r°7r°7 
and 7r 7r°7T where, respectively, one and two photons are 
not observed). The background component of the fitted 
function is shown as a gray surface. The background is 
estimated and subtracted as a function of z as indicated 
in Figure 7 and amounts to (2 — 2.5)% for z < 0.3 and 



(1.5 - 2)% for 0.3 < z < 0.6. 

A similar procedure is used in Figure 5 for determining 
the background at the ir°r) peak, which amounts to a total 
of 11% over the range —0.6 < z < 0.6 and is distributed 
as shown in Figure 7. 

Figure 6 shows the 7777 peak and its fitted background, 
amounting to a total of 8% over the range < z < 0.6 
with the angular distribution shown in Figure 7. 

IV. THE CROSS SECTION 

A paper on the 7r°7r° reaction has been published [2]. 
The measured pp — > 7T°7r° differential cross section is 
shown in Figure 8 for 3 of the 17 energy points of the data 
sample and over an angular range limited to < z < 0.6 
by the detector acceptance. The cross section integrated 
to various values of z max is shown in Figure 9 for all 
energy points. 

Figures 8 and 9 also show a binned maximum- 
likelihood fit to the cross section. The fit is performed 
simultaneously on all 17 energy points and over < 
z < 0.6. Within this range, the number of background- 
subtracted 7T°7T° events is 431,625. The parameteriza- 
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FIG. 4: The region of the peak (truncated at about 2% 

of its height) with a fit to the background shown as a gray 
surface. 




FIG. 5: The region of the peak (truncated at about 20% 
of its height) with a fit to the background shown as a gray 
surface. 

tion of Equations (6)-(8) is used setting J max = 4; the 
number of partial waves to include and their energy de- 
pendence were determined by searching for significant 
improvements of the x 2 - The mass and width of the 
Xco are constrained to the values (see Table IV) deter- 
mined by studying the process pp — > J/ip JlJ/ip — > e + e~ 
[6]. A magnified version of the plot of da/dz at E cm — 
3415 MeV, as well as a detailed description of the fit and 
its different components, is provided in [2] . The fit results 




FIG. 6: The region of the rjrj peak (truncated at about 12% 
of its height) with a fit to the background shown as a gray 
surface. 

are given in Table III. 

The fit presented in this section demonstrates the gen- 
eral structure of the pp — > ir°ir° angular distribution and 
estimates the number and amount of the contributing 
partial waves. However, this fit is not very sensitive to 
the value of the resonant amplitude Ar. The reason is 
that the size of the (interference-enhanced) resonant sig- 
nal is significant only at small values of z, while the fit 
is dominated by the high statistics of the nonresonant 
forward peak. The next section describes how the ex- 
traction of Ar is carried out. 

V. EXTRACTION OF 

B( X co -» PP) X B( X co 7r°7r°) 

As discussed in Section I (see in particular the dis- 
cussion of Equation 11), the natural region to exploit in 
order to obtain a model-insensitive measurement of the 
resonance amplitude Ar is the region at small values of z. 

To find the value of z at which the nonintcrfcring con- 
tinuum starts to play a significant role, we perform inde- 
pendent fits in each bin Az at small z. Each fit uses the 
parameterization of Equation (6) with Be l&B set equal to 
zero. The parameterization used for the helicity-0 com- 
ponent of the continuum is 

A 2 = a + a\x + a 2 x 2 , (12) 

to reproduce the moderate energy dependence of the 
cross section, see Figure 9. 

Since the growth with z of the helicity-1 component is 
not accounted for in these fits, the estimate of Ar, which 
in reality is a constant, will show an artificial decrease at 




FIG. 7: The number of candidate events (Tot) as a function of z at the XcO mass, the estimated background events (Bkgd), and 
a polynomial fit (solid curve) to the number of background events. Corrections for detector acceptance and efficiency have not 
yet been applied. For the r/rj channel, the background distribution was independent of energy and the data from all 17 energy 
points have been merged to reduce bin-to-bin fluctuations. 




FIG. 8: The pp — > -k°-k differential cross section versus z = | cos#*| at different E cm . A fit using Equations (6)-(8), with values 
from Table III, is shown. 



values of z where the helicity-1 component can no longer 
be neglected. The value of Ar can then be derived by 
using all the data up to some maximum value of z identi- 
fied by the fits to the individual bins. The fit results for 
Ar from individual bins of Az up to z = 0.3 are shown 
in Figure 10-top and the values of An derived by using 
the data within < z < z max are shown in Figure 10- 
bottom. 

The observed drop in An of Figure 10-top at z > 0.15 
shows that the helicity-1 component is no longer negli- 
gible at this value of z. We therefore restrict our region 
to 

z < z max = 0.125 (13) 
The restricted z range increases the statistical uncer- 



tainty on Aji, but ensures there is little effect on An from 
systematic uncertainties in our knowledge of B. The fit is 
described in detail in [2] and is reproduced in Figure 11. 

The value obtained for the resonant amplitude An, is 
given by: 

A\ = ttA 2 x B(xcO -> PP) x B( Xc0 - ttV) , (14) 

where A is the center-of-mass de Broglie wavelength, 
which gives: 

B(XcO -> pp)xB{ X cO -> ttV) = (5.09±0.81±0.25)xl0~ 7 . 

(15) 

The uncertainties are statistical and systematic, respec- 
tively. The dominant systematic errors arise from the 
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TABLE III: Fit results for coefficients and phases of the partial- wave expansion of Equations (6)-(8) for the 7r°7r°, -K°rj and r\r) 
channels. A linear energy dependence is included when found necessary. The errors are statistical. 
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luminosity determination (2.5%) and the knowledge of 
M Xc0 from [6] (2.5%). The systematic error due to the 
uncertainty in the helicity-1 continuum \B e l&B \ 2 is 1.2%. 

The phase between the helicity-0 nonresonant ampli- 
tude A and the resonant amplitude Ar is 6a = (39 ± 5 ± 
6) degrees. In the region < z < 0.125 of the fit, no 
dependence of 6 a on z and x is found. The values of So, 
62, and S4 are given in Table III. 



VI. THE -qr] CROSS SECTION 

The measured pp — > 7777 differential cross section is 
shown in Figure 12. Figure 13 shows the integrated cross 
section versus E cm . As in the 7r°7r° case, the cross section 
is dominated by the nonresonant continuum pp — > 7777, 
which has a smooth dependence on the energy The 7777 
channel also shows a resonance signal near the mass of 
the XcO- The interference pattern is different from the 
one observed in the 7r°7r° channel. There is destructive 
interference on the low-energy side of the resonance and 
constructive on the high-energy side, with the resonance 
peak shifted to above the XcO mass. 

As in the 7r°7r° analysis, a binned maximum-likelihood 
fit to the pp — > 7777 differential cross section is performed 
simultaneously on all 17 energy points and over an angu- 
lar range limited to < z < 0.6 by the detector accep- 
tance. Within this range, there are 19,675 background- 
subtracted 7777 events. The parameterization used is given 
in Equation 6, with the partial- wave expansion of Equa- 
tions 7 and 8. The mass and width of the XcO ar e con- 
strained to the values determined by means of the J/ip 7 
channel [6]. 

The result of the fit is shown in Figures 12 and 13. 
As in the 7r°7r° case, the contribution of the "pure" res- 
onance, A 2 R /{x 2 + 1), is negligible. The B 2 (helicity-1) 
term dominates at larger z, but its growth begins at a 
larger value of z than in the 7r°7r° case. The interfering 



helicity-0 continuum, of magnitude A 2 , dominates for a 
large part of the angular range (0 < z < 0.5) and provides 
the amplification for the interference-enhanced resonance 
signal seen in the separation between da/dz and A 2 +B 2 . 
Since the number of events available in the 7777 channel is 
limited (about 20 times fewer than in the 7r°7r° channel), 
the angular distribution within the available range in z 
is not as well resolved as in the 7r°7r° and ir°i~i channels. 
Compared to the 7r°7r° and 7r°?7, the 7777 channel requires 
a smaller number of parameters. In the best fit (see Ta- 
ble III) Co and C2 have a linear energy dependence; C4 
vanishes and is subsequently constrained to zero; and 
C21 and C41 do not require an energy dependence. The 
phases So, S2 and 64 are equal to each other, different 
from zero (that is the phase of the resonant amplitude is 
different from the phases of the nonresonant amplitudes) 
and constant in energy. The phases of the helicity-1 con- 
tinuum (only the difference between Si and 8\ is measur- 
able) also are not significantly different from each other 
and do not exhibit an energy dependence. Note that even 
in the 7r°7r° and 7T°77 cases the phases of the nonresonant 
amplitudes of a given helicity are very similar to each 
other. 



VII. EXTRACTION OF 

B(xco — PP) X B( X co — 1777) 

As in the tt°tt° analysis, the effect of the helicity-1 com- 
ponent of the continuum is investigated by performing a 
series of independent fits on each bin Az with Be lSB set 
equal to zero. Figure 14-top shows the fit results for Ar 
in every bin Az up to z = 0.4 . No drop of Ar is ap- 
parent in the region < z < 0.4 , indicating that any 
significant rise of B 2 occurs at larger z than in the 7r°7r° 
case. The fit in Figure 12 also shows that B 2 is small 
compared to A 2 below z w 0.4 . This is helpful because, 
in order to statistically resolve a clear resonant signal in 
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FIG. 9: The measured pp 
over < z < plotted versus E, 

tions (6)- (8) is shown. 



7r°7r° cross section integrated 



A fit using Equa- 



the plot of the 7777 cross section versus E cm , it is neces- 
sary to integrate over a larger z-range than in the 7r°7r° 
analysis. Using a larger z-range produces a slightly larger 
uncertainty in the relative contributions of A 2 and B 2 , 
and consequently in Ar. 

As in the 7r°7r° analysis, the information from the dif- 
ferent bins can be merged by performing new fits, inte- 
grating over ranges < z < z max with increasing z max . 
Equation 6 is used. The helicity-1 continuum is con- 
strained to the estimate from the partial wave expansion 
fit of Figure 12. The parameterization of the magnitude 
of the helicity-0 continuum is 



0.15 0.2 

z 



FIG. 10: Fits for the resonant amplitude of pp — > XcO — > 7r°7r°. 
Top: the fits are to the data of each bin in Az — 0.025. B 
is set to zero (see text). Bottom: each fit is performed over 
a range < z < z max for increasing z max and \B e tSn \ 2 is 
inserted from the fit of Figure 8. (The inner error bars are 
statistical while the outer ones indicate the uncertainty on 
C 2 \ C\ and <5l-<5^ in Table III). 



where the necessary parameters are added, with increas- 
ing z max , by searching for significant improvements in 
the x 2 °f the fit. The fact that an energy/polar-angle 
mixing term {asz 2 x) is required is already noticeable in 
Figure 13. The phase 8a does not require any significant 
dependence on x or z. The results for the resonant ampli- 
tude An as a function of z max from these fits are shown 
in the bottom of Figure 14. The estimate of Ar is stable 
for z max > 0.2 and the value of Ar is extracted from the 
fit performed over the angular range < z < 0.35 . The 
result of this fit is shown in Figure 15. The angular region 
of integration is larger than in the tt°tt° analysis, but still 
guarantees that the nonintcrfcring helicity-1 continuum 
(B 2 ) is relatively small. 

An equation analogous to Equation 14 gives 



:° : 5 3 ) x 10- 7 



A 2 = ao + a\x - 



2 1 2 1 4 
- a,2Z + C13Z x + CI4Z , 



(16) 



B(Xc0 -> PP) x B( X c0 - W) = (4.0 ± 1.2 

(17) 

The uncertainties are statistical and systematic, respec- 
tively. The dominant systematic error arises from the un- 
certainty in the helicity-1 continuum \B e i&B \ 2 : {tl'_t%)- 
The phase difference between the helicity-0 nonreso- 
nant amplitude A and the resonant amplitude Ar in the 
angular region < z < 0.35 is 5a = (144±8±6) degree; 
no dependence on x and z is observed. The phase 5a, 
responsible for the shape of the interference pattern, is 
different from the corresponding phase in the 7r°7r° chan- 
nel. In the 7777 channel, the interference is destructive on 
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FIG. 11: The pp — > TT°ir° cross section integrated over < 
z < 0.125 plotted versus E cm . The fit of Equation 6 is shown 
along with its components. The contribution of \B e lSB \ 2 is 
constrained to the estimate obtained from the partial-wave 
expansion fit of Figure 8, and its integration corresponds to 
the separation between the curves J (A 2 + B 2 ) dz and J A 2 dz. 
The "pure" Breit-Wigner curve, J[A%/(x 2 + 1)] dz, is shown 
multiplied by a factor 20. 



the low-energy side of the resonance and constructive on 
the high-energy side. 



VIII. ESTIMATE OF M xco AND T xco USING 
THE 7T 7r° AND 7777 SAMPLES ALONE 

The results presented so far have been obtained with 
the Xco mass (M Xe0 ) and width (T Xc0 ) constrained to 
the high-precision values coming from the analysis of the 
J/tp 7 channel [6] . It is worth checking the capability of 
the interference technique in determining M Xc0 and r Xc0 , 
in addition to determining the product of the input and 
output branching ratios. The fits shown in Figures 11 and 
15 are redone in Figure 16, with M XcQ and r Xc0 left as free 
parameters. The results for M Xc0 , T Xa0 , B in x B out and 
5a for both 7r°7r° and 7777 channels are given in Table IV 
in Section XI. 

As can be seen, the interferometric technique is able to 
determine the resonance parameters without relying on 
additional inputs. The uncertainties on M Xc0 and r Xc0 
from the analysis of the 7r°7r° channel alone are larger but 
still comparable to those from the analysis of the J/ipj 
channel. 



IX. THE 7r°77 CROSS SECTION 

The measured pp — > 7r°7i differential cross section is 
shown in Figure 17 for 3 of the 17 energy points over the 
angular range —0.6 < z = cos#* < 0.6. The integrated 
cross section versus E cm is shown in Figure 18; there are 
85,751 background-subtracted ir°r) events. Also shown 
is a binned maximum-likelihood fit performed simulta- 
neously on all 17 energy points to Equation 6-8. The 
parameter Ar is set to zero, as the decay of any charmo- 
nium state into w ri is suppressed by isospin conservation. 

Since the initial pp state is a linear combination of C- 
eigenstates, the final 7r°7i angular distribution must be 
symmetric in the center-of-mass polar angle. Figure 7 
shows that the data, uncorrected for acceptance, are ob- 
viously not symmetric for \z\ > 0.5. This arises because 
slow moving 77's, i.e. those emitted backward in the cen- 
ter of mass, have a significantly larger opening angle for 
the decay photons than 7r°'s of the same center-of-mass 
angle and thus a greater probability of producing a pho- 
ton outside the acceptance, particularly near the accep- 
tance boundaries. The calculated acceptance is shown in 
Figure 3. Figure 17 shows that the forward/backward 
symmetry is recovered in the acceptance-corrected angu- 
lar distribution. 

The values of the fit parameters arc given in Tabic III. 
The need for the J = 4 amplitude is more evident than 
in the 7r°7r° channel and a fit with J max = 2 cannot re- 
produce the multiple minima at z — and \z\ ~ 0.45. 
The dip in the differential cross section at z = is 
very pronounced; see Figure 17. As explained before, 
the hclicity-1 component B 2 vanishes at z — 0, and the 
hclicity-0 component A 2 of the cross section is very small 
at z = due to a local cancellation of the involved par- 
tial waves. The pronounced minimum in the differential 
cross section is not present everywhere within the E cm 
range 2911 MeV to 3686 MeV [7]. 

Fits to the cross section are performed with the reso- 
nance amplitude Ar as a free parameter. The procedure 
to obtain the upper limit 

B(Xc0 - Pp) x B( X c0 -> A) < 4 x 10- 8 (18) 

at 90% confidence level is described in detail in [3] . This 
upper limit is one tenth of the values for 7r°7r° and 7777 
given in Table V. 

X. THE ttV AND 7777' CROSS SECTIONS 

The small 77' — > 77 branching ratio and a larger back- 
ground limit the achievable precision of the study of chan- 
nels (4) and (5). The presence of 7r r/ and 777/ events can 
be recognized in Figure 2. 

The event selection and variables employed are sim- 
ilar to those described in Section II. Events with four 
CCAL energy deposits greater than 100 MeV are selected 
and a 5% confidence level on a 4C fit to pp — ► 47 is re- 
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FIG. 12: The pp — > rjrj differential cross section versus z = | cos#* at different E cm . To reduce statistical fluctuations, the 17 
energies are displayed merged as 3339.5 MeV < E cm < 3400.1 MeV (left), 3406.1 MeV < E cm < 3430.1 MeV (center) and 
E cm = 3469.9 MeV (right). A fit using Equations (6)-(8) is shown along with its components (see Table III for values of fit 
parameters) . 



quired. The event topology is denned as the combina- 
tion of the four photons [named as 7r° (or 77) — > 7172 and 
i]' ~^ 7374] w ith the highest confidence level of a 6C fit 
to pp — > ir°r)' (or 7777') — * 47. Coplanarity and colinearity 
cuts as listed in Table II are then applied. For both anal- 
yses, it is additionally required that m li7j > 250 MeV 
for combinations i,j = 1,3; 1,4; 2,3; and 2,4 (i.e. 
the photon pairings not chosen as the event topology) 
to reject ir° contamination. For the 777/ analysis, it 
is additionally required that the sum of the "wrong" 
paired combinations m 7l73 +m 7274 and 
greater than 2.5 GeV. This cut does not seriously affect 
the acceptance as seen by MC simulation. A substan- 
tial improvement of the resolution of the 77' peak (from 
(7 — 40 MeV to 16 MeV) is obtained by using the output 
values of the photon energies and positions of a 5C fit to 
pp — ► 7T° (or 77) 7374 — ► 7172 7374. The resulting distri- 
bution of shown in Figure 19 for 7r°?7' and 7777' 
events. A cut \m 73l4 — m,y \ < 40 MeV is applied. 

To determine the angular distributions with adequate 
resolution, all 17 energy points are merged together. A 
plot of m 7l72 versus 7n 7374 (prior to the 5C fit) shows 
that the "berm" along the line m 7374 =m^' is negligi- 
ble compared to the one along and along 
Hence, the background is simply determined 
by fitting 1-dimensional projections (after the 5C fit) 
like those in Figure 19. This is done for each bin Az, 
where z = cos 9* (z = cos 0*) for the 7r°?7' (7777') anal- 
ysis. In the range \z\ < 0.6 there are 15,097 candidate 
7t t/ events of which ~ 33% are estimated to be back- 
ground. In the range \z\ < 0.6 there are 1166 7777' can- 
didates with an estimated background of ~ 25%. As 
in the 7r°77 case, the acceptances for n°r]' and 7777' are 
not forward/backward symmetric. The symmetries are 
recovered in the background-subtracted and efficiency- 
corrected angular distributions as shown in Figure 20. 



Figure 21 shows the pp —> ir°r]' and pp — > 7777' cross 
sections integrated for \z\ < 0.3 versus E cm (the back- 
ground has not been subtracted in the plot). In prac- 
tice, the background has been determined at each energy 
point with fits like those in Figure 19. The data are not 
sufficient to perform partial-wave analyses to distinguish 
the hclicity-0 and helicity-1 components of the continuum 
and the fits in figure 21 are performed simply on the in- 
tegrated cross section. Equation 11 is used. Given the 
small values of \z\ in these fits, the helicity-1 component 
B 2 is fixed to zero for both channels. The background 
is parameterized as a polynomial and fit simultaneously 
with the total cross section. 

Charmonium decay into 7r°7/ (7777') violates (satisfies) 
isospin conservation. For the ir rj' channel, the quality of 
the fit remains unaltered by allowing the resonant ampli- 
tude A free or fixing it to zero. The x 2 is 23.9 in both 
cases and the degrees of freedom are 27 or 29, respec- 
tively. The 90% confidence level upper limit is 

B(Xc - PP) x B( X c0 - A') < 2.5 x lO" 7 . (19) 

For the 7777' channel, the fit to the total cross section 
changes from the dashed to the solid line in Figure 21 
when a resonant amplitude A is allowed. The \ 2 /DOF 
goes from 30.4/31 to 25.6/29. If the feature observed in 
the data is due to interference between the XcO and the 
continuum, the solid line fit would imply 

B( X cO -> PP) x B(xcO -> W) - (2.1± 2 ; 3 5 ) x 10- 6 . (20) 

The error quoted is statistical and dominates the system- 
atic uncertainty. 
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FIG. 14: Fits for the resonant amplitude of pp — > XcO —* 
rp). Top: each fit is performed in a bin of size Az = 0.1 
independently from the other bins; Equation 6 is used with 
the helicity-1 continuum \B e lSB | 2 fixed to zero. Bottom: 
each fit is performed over a range < z < z m ax for increasing 
Zmax, Equations 6 and 16 are used, while \B e lSB \ 2 is taken 
from the partial-wave expansion fit of Figure 12 (the outer 
error bars show the uncertainty due to the errors on C\ and 
Cl in Table III). 



For 7r°7r° and 7777, the phase 5 a at small z between the 
hclicity-0 nonresonant and resonant amplitudes is also 
given. There is good agreement for the values of M Xc0 
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FIG. 13: The measured pp — > r/rj cross section integrated over 
< z < Zmax plotted versus E cm - A fit using Equations (6)- 
(8) is shown. 



XI. RESULTS 
A. The Xco State 

The E835 measurements of the XcO resonance are sum- 
marized in Tables IV and V. The results from the chan- 
nels J/ipj and n°TT have been already published in [6] 
and [2], respectively. 

The header section of Table IV indicates the common 
formation channel (pp) and the different \c0 decay chan- 
nels: J/ip'y, n°TT°, and 7777. The bottom part of the table 



XcO ' x XcO ' 



XcO 

sample of 7777 events gives lower-precision measurements; 
the width is in agreement with the other two channels, 
while the mass is slightly underestimated (correlated with 
a probable slight overestimate of the phase 5a)- The 
background from misidcntificd events from different pro- 
cesses is negligible for all three cases. There is an ab- 
sence of nonresonant production in the J/ip 7 channel. 
The large nonresonant continuum in the 7r°7r° and 7777 
channels, although useful as the provider of the ampli- 
fication for the interference pattern, requires additional 
parameters to be fit. In particular, the phase 5a and the 
size of the interfering continuum increase the coupling 
among the fit parameters. As a result, the statistical 
uncertainties are larger than in the J/ip 7 channel. 

Table V presents the results for Bi n x B out and the 
phase 5a as determined for the neutral pseudoscalar 
channels by constraining M Xc0 and r Xc0 to the measure- 
ments from the J/ip 7 channel. The J/ip 7 data sample 
was collected at the same time as the four photon events. 
The energy and luminosity determinations are the same 
and Bi n x B out as for all these channels. By using the values of M Xc0 and 

out 



reproduces the results for M- 

determined by each of the three analyses independently. r Xc0 from the J/ip'y analysis, the values of Bi n x B, 
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TABLE IV: E835 Results for the XcO with M Xc0 , T Xa0 , and B in x B out as free parameters. 



B in = 
B ou t = 


Commo 

B( X co^ J/i>i) 


n channel B(\ c o 
S( X cO - ttV) 


->pp) 

B(XcO -> w) 


M Xc0 (MeV/c^) 
Txco (MeV) 
S in x Bout (1(T 7 ) 
5a (degree) 


3415.4 ±0.4 ±0.2 
9.8 ±1.0 ±0.1 
27.2 ± 1.9 ± 1.3 


3414.7^'i ± 0.2 

8.6ti 3 ±0.1 
5.42±g"|^ ± 0.22 
47 ± 10 ± 6 


3412. 2t'ii ± - 2 
10.3+f-? ±0.1 

+ 1.2 +0.5 
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173t\l ± 6 
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FIG. 15: The pp — > r/rj cross section integrated over < 
z < 0.35 plotted The fit of Equation 6 and 16 

is shown along with its components. The \B e l * B | 2 term is 
constrained to the estimate from the partial-wave expansion 
fit in Figure 12; its contribution is shown by the separation 
between the curves / (A 2 + B 2 ) dz and / A 2 dz. The "pure" 
Breit-Wigner curve, f A\j (x 2 ± 1) dz, is shown multiplied by 
a factor 10 to make it comparable to the observed interference 
signal. 



and 8a in the 7r°7r° channel are determined with higher 
accuracy and precision. The values in Table V arc our 
final results for B in x B out and 5 a- 

The value of the phase 5a between the helicity-0 non- 
resonant and resonant amplitudes is determined in a re- 
stricted angular region at small z (z < 0.125 for 7r°7r°, 
z < 0.35 for r\r\ and z < 0.3 for rjij'). No appreciable en- 
ergy dependence is seen within the considered regions of 
z. The phase 5a is responsible for the shape of the inter- 
ference pattern seen in the cross section and is reasonably 
well measured. Its specific value is determined by the lo- 
cal combination of several partial waves (see Equation 7), 
which largely cancel each other. 

In the ratio between Bi n x B out of two analyzed chan- 
nels, the common Bi n and some systematics cancel out. 
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FIG. 16: The ir°7T° (left) and r\r] (right) cross sections inte- 
grated over the indicated angular ranges and plotted versus 
E cm . The fits of Figures 11 and 15 are re-performed allowing 
M Xc0 and r Xc0 to be free parameters. 



We obtain 

Bjxco -> ryq) _ 27+0 .io m ] 

B(Xc0 — y Tr®iT®) 0.79 ± U.27_q Qg (21) 

The BES experiment has reported the measurement of 
B(Xc0 -> 7r°7r°) = (2.79 ± 0.32 ± 0.57) x 10~ 3 , which 
agrees with the 1985 measurement from Crystal Ball [1]. 
However, a consistent and more precise determination of 
this quantity may be computed by using isospin sym- 
mctry: B( Xc0 -> vT°7r°) = \B{ Xc0 -> tt+ij-) = (2.5 ± 
0.35) x 10~ 3 , where B( Xc o tt+tt-) = (5.0 ±0.7) x lO" 3 
is taken from [1]. This and Equation 21 provide 

B(xco -» rm) = (1.98±0.68 stat t°-^ s ^±0.28 PI5G )xl0- 3 , 

(22) 

where the subscript PDG labels the uncertainty derived 
from errors listed in [1]. This is in agreement with a 
measurement B( Xc a -> w) = (2.02 ± 0.84 ± 0.59) x 10" 3 
reported by BES [8], again in agreement with a 1985 
measurement from Crystal Ball [1]. 



Non-resonant pp Annihilation into Two 
Pseudoscalar Mesons 



In addition to the charmonium results, this work pro- 
vides fits to the cross sections for three pseudoscalar- 
pseudoscalar meson states in antiproton-proton annihi- 
lations in the energy range 3340 MeV —3470 MeV (see 
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FIG. 17: The pp — > iv r) differential cross section versus z = cosd^ at different E cm . A fit using Equations (6)-(8) with Ar = 
is shown along with its components. The values of the fit parameters are reported in Table III. 



TABLE V: E835 Results for the XcO with M Xc0 and r Xc0 constrained to the values of the J/V>7 results in Table IV. 



-Bin — 
B ou t = 


B( X cO - ttV) 


Commo 
B(XcO -> rm) 


q channel B(\ 
B(XcO -> W) 




B(XcO -» A') 


Bi n x B out (10" Y ) 
8a (degree) 


5.09 ±0.81 ±0.25 
39 ± 5 ± 6 


4.0 ± 1.2+;;^ 

144 ± 8 ± 6 


zl -15 

189 ± 20 


< 0.4 (90%CL) 


< 2.5 (90%CL) 



Table III) . The differential cross section provides insights 
on the dynamics of these processes at the examined en- 
ergies. The differential cross sections as a function of 
the production angle of the meson pair are shown in Fig- 
ures 8, 17 and 12 for tt°tt°, ir r/ and 7777, respectively, and 
in Figure 20 for tt°t]' and 7777'. Tables with the numerical 
values for tt°tt°, Tr°r] and r\r\ can be found in [3]. 

The partial-wave expansion fits described in Sec- 
tions IV (7r°7r°), IX (w°r]) and VI (7777) indicate that 
the size of the angular-momentum contribution decreases 
rapidly with J and only J = 0, 2 and 4 are necessary to 
describe the angular distributions. The statistical reso- 
lution for 7r ?/ and 7777' is limited. However, their angular 
distributions are well fitted by an even-power expansion 
up to z 8 , hence are consistent with J ma x = 4. 

Applying the semi-classical relation hLp p ~ bpi 
(where L p p is the orbital angular momentum of the pp 
system, b is the impact parameter, and pi is the initial 
state center-of-mass momentum) , it can be inferred that 
a larger number of partial waves would participate. This 
is probably true for other reactions and is known to be 
true for elastic pp scattering. However, the pseudoscalar- 
pseudoscalar meson states selected here are unlikely to 
be produced in peripheral impacts, since they require the 
annihilation of one, two or three valence qq pairs of the pp 
initial state. If only one qq pair annihilates, the remain- 
ing valence quarks of the proton (qq) and antiproton (qq) 
must separate and rearrange themselves into qq + qq . 
This is unlikely to happen at large impact parameters, 
where the partons tend to retain their original large longi- 
tudinal components of momentum. An indirect observa- 



tion of a sharply decreasing trend in the relative contribu- 
tion of increasing L pp is provided by the exclusive process 
pp — ► cc, which necessarily requires total valence-quark 
annihilation. The ratios Tp p /T g \ uons for the charmonium 
states that couple to L Pp = (rj c , J/ip and ip') are 5-10 
times larger than for those that couple to L pp = 1 (the 
Xcj's) [1]. Table I indicates that only odd values of Lp p 
can produce the J PC = cvcn ++ of a pseudoscalar me- 
son pair, and each L Pp feeds both J = L Pp ± 1 . It is 
then reasonable to expect larger contributions from the 
partial waves with J — and 2 as compared to J = 4 
as observed in Table III. It may also be noticed that the 
estimated differences of phase between amplitudes with 
the same helicity are relatively small. 



XII. DISCUSSION AND CONCLUSIONS 

E835 has studied the formation of the \c0 state of char- 
monium in antiproton-proton annihilation and its subse- 
quent decay into pseudoscalar-pseudoscalar mesons. In 
the 7r°7r° and 7777 channels, an interference-enhanced pat- 
tern is evident in a cross section dominated by the non- 
resonant production of pairs of pseudoscalar mesons. 

The choice of performing this study on the XcO reso- 
nance is a consequence of the J PC = ++ quantum 
numbers of the \ c o, which allows it to decay into a pseu- 
doscalar meson pair. The primary goal of E835 was the 
determination of the resonance parameters through the 
study of the process pp — > J/tp 7, J/ip — > e + e~ to com- 
plete the program of studying the %cj triplet initiated 
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FIG. 18: The measured pp — > tv r\ cross section integrated 
over \z = cos 6„ \ < z max plotted versus E cm . A fit using 
Equations (6)-(8), with Ar = 0, is shown. 



by E760. Thanks to the antiproton source developments 
mentioned in the introduction, a large integratged lumi- 
nosity was collected in the Xco region in the year 2000 run 
of E835. The presence of a large nonresonant continuum 
cross section for the final 7r°7r° and r]T] states, as compared 
to the resonant production through a charmonium inter- 
mediary state, was known at the outset. However, the 
awareness that the interference mechanism would pro- 
duce an enhanced interference pattern in the cross sec- 
tion, hopefully large enough to be detected, was a strong 
motivation to pursue the analysis. Of course, the physics 
behind the interference mechanism was well known prior 
to this work. What is innovative is the exploitation of 
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FIG. 19: The r)' (958) 7374 peak for the tt ?/ (left) and 
r/rj' (right) selections. The peak at lower energy is due to 
cj(782) — > 7r°7, 7r° — > 77 events where one of the 7r° decay 
photons is not observed. The fit is a gaussian plus a first- 
degree polynomial. Arrows indicate the applied cut. 
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FIG. 20: The pp — > 7r°r/' and pp — > 7777' differential cross sec- 
tions at E cm — 3340 — 3470 MeV. The background has been 
subtracted. A fit to a power expansion in even-powers up 
to z 8 (corresponding to J max — 4) is shown to verify for- 
ward/backward symmetry. 



such a mechanism to detect and measure a resonant sig- 
nal that would, if the interference could be turned off, be 
almost two orders of magnitude smaller than the nonres- 
onant cross section. Several other measurements of in- 
terference patterns exist in particle physics and many in 
nuclear physics, but the usual experience is with a larger 
resonance amplitude and a smaller interfering continuum. 
This specific analysis requires the correct separation of 
the two different continuum components, one interfering 
and the other not interfering with the resonance. The 
effectiveness of this analysis has been demonstrated. 

A reason for pursuing this study was the search for 
alternative means of discovering and measuring charmo- 
nium states and possible hadromolecular states. Now 
that confidence has been gained that measurements 
of resonances can be accomplished in hadronic decay 
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channels where the nonresonant production of the fi- 
nal state dominates, new strategies can be considered. 
For example, other than the above mentioned search for 
hadromolecules, the study of the poorly known charmo- 
nium singlet states could be performed by investigating 
hadronic final states, relying on the enhancement pro- 
vided by the interference. 
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FIG. 21: The pp — ► ir°i~i' and pp — > r\rj cross sections inte- 
grated over \z\ < 0.3 (z = cos#* and z = cos 91 for tt rf and 
r\r) , respectively) plotted versus E cm . In these two plots, the 
amount of cross section due to background (Bkgd) has not 
been subtracted from the total (Tot), and the fits shown are 
performed simultaneously on Tot and Bkgd cross sections (see 
text). 
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